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The formation of the coronal flow/ADAF 
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Abstract. We develop a new method to describe the accretion flow in the corona above a thin disk around a black 
hole in vertical and radial extent. The model is based on the same physics as the earlier one-zone model, but 
now modified including inflow and outflow of mass, energy and angular momentum from and towards neighboring 
zones. We determine the radially extended coronal flow for different mass flow rates in the cool disk resulting 
in the truncation of the thin disk at different distance from the black hole. Our computations show how the 
accretion flow gradually changes to a pure vertically extended coronal or advection-dominated accretion flow 
(ADAF). Different regimes of solutions are discussed. For some cases wind loss causes an essential reduction of 
the mass flow. 
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1. Introduction 

The new fascinating results from XMM-Newton and the 
Chandra X-ray Observatory allow a deeper insight into 
the physical processes in many astrophysical objects, on 
scales ranging from binary stars to galaxies. One of these 
topics is the accretion onto black holes. This can happen 
in form of an advection-dominated accretion flow (ADAF) 
or via a standard geometrically thin Shakura-Sunyaev ac- 
cretion disk. Since an ADAF is only possible in the inner 
region around the central accretor one has in many ob- 
jects an advection-dominated accretion flow in the inner 
region, fed by the mass flow through a standard disk at 
larger distances from the black hole (for a recent review 
on ADAF models see Narayan 2002). The situation is the 
same for accretion onto black holes in galactic X-ray bina- 
ries and onto supermassive black holes in active galactic 
nuclei (AGN). 

The spectra arising from these two different ways of 
accretion are very different, a soft multi-temperature spec- 
trum from the black-body like accretion disk and a hard 
power law spectrum from the Comptonizing much hotter 
gas in an ADAF (one observes the combination of both 
originating at different distances from the black hole). The 
observed change of the spectral type, the soft/hard spec- 
tral transitions in X-ray binaries, were successfully mod- 
eled as related to a changing mass flow rate by Esin et al. 
(1997, 1998), thereby supporting this picture of an ADAF 
in the inner region surrounded by a standard disk. This 
means that at a certain distance from the black hole the 
accretion mode changes from the flow via a thin disk to 
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a hot gas flow. The distance where this happens depends 
on the mass flow rate from outside. 

Meyer, Liu and Meyer-Hofmeister (2000) developed a 
model for a corona above a geometrically thin standard 
disk around a black hole (basically the same physics as al- 
ready discussed earlier for disks of dwarf nova systems, 
where the compact object is a white dwarf (Meyer & 
Meyer-Hofmeister 1994)). The corona is fed by gas which 
evaporates from the cool thin disk underneath. An equi- 
librium establishes between the cool accretion stream and 
the hot flow. The efflciency of this process increases to- 
wards the black hole. This means that at a certain distance 
all matter is evaporated and the disk is truncated. From 
this distance on inward all gas is in the hot flow and pro- 
ceeds towards the black hole as an advection-dominated 
flow. This original model is a simplifying description of 
the two-dimensional accretion flow (of vertical and radial 
extent with a free boundary condition on a radially ex- 
tended surface). The resulting evaporation efflciency as a 
function of distance (and of the central black hole mass) 
allows to determine the location at which the thin disk is 
truncated. 

For some X-ray binaries in quiescence the location of 
the inner edge of the disk can be deduced from the or- 
bital velocity there, which is inferred from the observed 
Hq line proflles. This requires that the disk temperature 
at the innere edge does not exceed the ionization tem- 
perature of hydrogen. Using the ADAF model for the fit 
to the spectra of soft X-ray transients as e.g. A620-00, 
V404 Cyg and Nova Mus 1991 the mass accretion rates 
were derived (Narayan et al. 1996, 1997). If one compares 
these results with the disk truncation radius from the 
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evaporation model the agreement is reasonable (Liu et al. 
1999, Meyer-Hofmeister & Meyer 1999). The change from 
disk accretion to an ADAF was also investigated for sev- 
eral low luminosity AGN and elliptical galaxies (Quataert 
et al. 1999, Gammie et al. 1999, Di Matteo et al. 1999, 
2000), also the theoretically expected disk truncation (Liu 
& Meyer-Hofmeister 2001). The results confirm that the 
truncation of the thin disk is located at smaller distance 
from the black hole for higher mass flow rates in the thin 
disk. Only for low luminosity AGN the observed spectra 
seem to demand a disk truncation at radii too small for 
the accretion rate (Quataert et al. 1999). But this appar- 
ent discrepancy might be resolved by the effect of magnetic 
fields from a dynamo in the underlying disk on the coronal 
gas flow (Meyer & Meyer-Hofmeister 2002). 

We now present work which is an essential step be- 
yond the one-zone model. We develop a new method to 
describe the accretion flow in the corona above a thin disk 
in its vertical and radial extent. The model is based on the 
same physics as in the one-zone model, but now modified 
including inflow and outflow of mass, energy and angular 
momentum from and towards neighboring zones. In the 
earlier model inflow from outward regions was neglected. 
But this is necessary if one considers regions inside the 
evaporation maximum (in the one-zone model at about 
300 Schwarzschild radii, compare Fig. 3 in Meyer et al. 
2000). For the inner regions it is also important to take 
into account different ion and electron temperature and, in 
the case of a high mass flow rate in the thin disk, the effect 
of Compton cooling of coronal electrons by disk photons. 
Both were included in the one-zone model description by 
Liu et al. (2002), where emphasis was put on evaluation 
of the coronal flow in the case of high accretion rates as 
in narrow-line Seyfert 1 galaxies. 

The co-cxistcncc of hot and cold gas around galactic 
black holes and in AGN was also investigated by Rozanska 
and Czerny (2000a, b). Different ion and electron temper- 
atures were already included. In their work the physical 
picture is basically the same as in ours, but the results 
differ in detail. The investigation focuses on the inner- 
most region near the black hole. In earlier work (Meyer- 
Hofmeister & Meyer 2001) we discussed the difference of 
the results. Processes of evaporation very close to the black 
hole were studied by Spruit & Deufel (2001). 

In Sect. 2 we discuss the physics of interaction be- 
tween the hot corona and a cool disk underneath. We in- 
troduce the modifications necessary for a consistent treat- 
ment of the mass, energy and angular momentum flow. 
We describe the procedure how one finds the solution 
for the two-dimensional accretion flow, put together from 
the separately computed structures for different radii, and 
we present the new computational results (Sect. 3). This 
multi-zone model describes the increase of the mass flow 
in the hot corona towards the black hole. A part of the 
matter is lost in a wind from the hot corona. In Sect. 4 we 
discuss the different regimes of solutions, wind loss from 
the coronal flow and the consequences of a varying mass 
flow rate in the thin disk on spectral transitions. 



2. The physics of the interaction of the hot 
corona and a cool disk underneath 

If we imagine hot gas above cool matter in a disk below 
both layers interact. This interaction happens in the grav- 
itational potential of the central star (as earlier considered 
near the white dwarf in binaries, a "siphon flow" , Meyer 
& Meyer-Hofmeister 1994). The hot corona conducts heat 
downward by electron conduction. At the bottom the tem- 
perature decreases from its high coronal value to a low 
chromospheric value, heat conduction becomes ineffective 
and the thermal heat flow has to be radiated away. The 
efficiency of radiation (in the optical thin case) depends on 
the square of the particle number density. If this density is 
too small the material will heat up and increase the den- 
sity in the corona. In this way the corona of a given tem- 
perature will "dig" itself so deep into the chromospheric 
layers that a density is reached which is able to radiate 
away the downward heat conduction and an equilibrium 
establishes between the cool accretion stream and the hot 
flow. The final density in the corona is determined by the 
pressure equilibrium at the interface. The hot flow contin- 
uously drains mass from the corona towards the central 
object. This is resupplied by evaporation from the surface 
of the cool disk as the corona tries to restore the density 
to the stationary level. 

2.1. The modeling of the corona above the thin disk in 
vertical and radial extent 

To describe the structure of the corona above the cool disk 
we take the standard equations of viscous hydrodynamics: 
conservation of mass, the equations of motion and the first 
law of thermodynamics. We want to determine the verti- 
cal structure of the corona at different distances r from 
the black hole with emphasis on the inner accretion re- 
gions. For the earlier "one-zone model" only the zone at 
the edge of the thin disk was considered. Since the evap- 
oration efficiency increases steeply in radial direction in- 
ward this approach gives already good results (but only 
for distances from the black hole where the evaporation 
efficiency indeed increases inward). We now want to in- 
vestigate the coronal structure in its full extension in two 
dimensions. We determine the vertically extended coro- 
nal structure for a series of successive radial zones. In this 
multi-zone model the divergence in radial direction is re- 
placed by inflow/outflow of mass and angular momentum 
in the zones (in the former one-zone model this could only 
be taken into account approximately). 

2.2. The equations 

We use cylindrical coordinates r, ip, z with the 2;-axis per- 
pendicular to the disk midplane. We consider stationary 
and azimuthally symmetric flows. We use basically the 
same equations as for the one-zone model, also used in 
Liu et al. (2002), but modified to take mass, energy and 
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angular momentum inflow and outflow in its radial depen- 
dence into account. 

In the following we list these five ordinary differential 
equations. The dependent variables are the vertical mass 
flow m = pVz {p density, vertical flow velocity), the ver- 
tical component of the heat flux Fc, the pressure P and the 
ion and electron temperature Ti and Te, the independent 
variables are r and z. 

We use the following equation of state 



P = mkTi + riekTe 



2po 



(1) 



with ni,n(.,Ti, Tg ion and electron number density and re- 
spectively temperature, pi = p is the ion mass density. 
We take a standard chemical composition (X = 0.75, Y ~ 
0.25) in the corona and the average molecular weight then 
po = 0.62. In a fully ionized gas of standard composition 
the electron number density is about 10% larger than the 
ion number density. The last expression for the total pres- 
sure in Eq. (1) is chosen for simplicity. For /xq = 0.62 it is 
exact if Te = Ti, but it overestimates P for Tg -C Tj by 
5%. We neglect this difference. 

Prom the conservation of mass we get our flrst equation 



2z 



r 



(2) 



with Vr the radial diffusive velocity and t^m the mass ad- 
vection modification term (see next section). The second 
term in Eq. (2) approximately takes into account the ef- 
fect of the changing channel cross section « (1 -|- z^/r^), 
for the ascending flow as its shape changes from cylindri- 
cal to spherical at large height z (compare Meyer et al. 
2000). 

The second equation is the z-component of the equa- 
tion of motion, 



pVz 



dvz 
dz 



dP 

dz 



GMz 



(r2 + z2)3/2 



(3) 



with G gravitational constant and M mass of the central 
black hole. The r-component of the equation of motion to 
first order gives the rotational velocity v^p as equal to the 
local Kepler velocity, = GMr^/{r^ + z^)"^/^, neglecting 
radial pressure gradients, small of order WTr/GMp < in 
our solutions. T stands for the temperature which deter- 
mines the hydrostatic layering, T ^ Ti. The (/^-component 
yields angular momentum conservation and determines Vr- 
The vertical conductive heat flux provides our third 
equation 



Fe = -KoT^/ 



■ dTe 
dz ' 



(4) 



with Ko = 10~^g cm s~^K~^/^ for a fully ionized plasma 
(Spitzer 1962). This formula for the heat conduction is 
valid if the electron mean free path is small compared 
to the length over which the temperature changes signifi- 
cantly. In the region of our solution where heat conduction 
is a significant contribution this is the case (occasionally 
at least marginally). 



The two remaining equations, the energy equations for 
ions and for electrons, as derived by Liu et al. (2002), in- 
chide the cooling and heating processes in the hot corona 
and the equations are now modified for the multi-zone 
model. For ions the energy balance is determined by vis- 
cous heating, cooling by collision with electrons and radial 
and vertical advection. The friction is taken proportional 
to the pressure with a standard a-prescription. 



-riPiVzUi) 
dz 



3 2 

-aPQ - Qie + m ■ -PiVrUi 

2 r 

2^ , . 



+ z^ 



(5) 



with the rotational frequency, pi ion density, a viscosity 
parameter. Ui is the specific energy of ions 



1^ 



7 



Pi. 



l-l Pt 



$), 



with $ the (Newtonian) potential 
GM 



$ = - 



(r2 + 02)1/2 • 



(6) 



(7) 



|aPri is the viscous heating rate per unit volume and g^e 
the rate of energy transfer from ions to electrons (Stepney 
1983) as taken in Liu et al. (2002). The energy advection 
modification term r]E is discussed in the next section. The 
contribution of the frictional stress to the divergence of 
the energy flow is already explicitly taken into account as 
frictional dissipation term in the equations. 

For the multi-zone model there are two changes com- 
pared to the one-zone model, (1) the mass and energy 
modiflcation factors and (2) we now take into account the 
radial dependence of in the formula for the speciflc 
energy (see Sect. 2.3.3). 

For electrons the energy balance is determined by 
the processes of heating by collisions with ions, cooling 
by bremsstrahlung (free-bound and free- free transitions), 
Compton cooling, and vertical thermal conduction. We ne- 
glect here the radial thermal conduction in a first approxi- 
mation. In the main part of the vertical structure solution 
it is not a dominant term (see Meyer et al. 2000). 

The second energy equation is 



d 
dz 



{peVzUe + -Fc) = lie " neniL{Te) - gcomp 



+m • -PeVrUe 

r 



2z 



j,2 _|_ ^. 



■{peVzUe + Fc), 



(8) 



where Ug is speciflc energy of electrons (defined corre- 
spondingly to Ui) and neniL{Te) is the bremsstrahlung 
cooling rate and gcomp the Compton cooling rate (Rybicki 
& Lightman 1979), as in Liu et al. (2002). 

In the energy balance for the electrons one could omit 
the terms multiplied with pe due to the low ratio of elec- 
tron to ion mass. We keep the equation for ions, Eq.(5), 
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and replace the equation for electrons by a joint equation 
for electrons and ions (with u the specific energy of elec- 
trons and ions). 

2.3. The modification terms for radial inflow/outflow 

We first derive the formula for the radial diffusive veloc- 
ity and then the mass and energy advection modification 
terms. For simplicity these modifications are taken in a 
height averaged way. 

2.3.1. The radial diffusive velocity 

The equation for local angular momentum conservation is 

^{2TTrpVrr^n) + ^{2Tirpv^r^Q.) + ^(27rrr™) = (9) 
ar oz or 

with the frictional stress Tr-,n 



= —aP = 



dfl 
dr 



3 

—-flu, 



(10) 



/X viscosity. For the lower, dominant part of the corona it 
is a reasonable approximation to neglect the z-dependcnce 
of fl (for an error estimate see Meyer et al. 2000). Mass 
conservation 



(11) 



(12) 



d d 
— {2TrrpVr) + T-{2TTrpv^) = 
or oz 

together with ^(r^fi) w yields 

d d 
2'nrpvr-^ir'^fl) + 2,-K—{r'^nu) = 0. 
oz or 

Using 

djr'^n) _ IrH} 
dr 2 r 

one finally obtains for the radial diffusive velocity (inward 
directed) 



^r = -^ 1 + 2^^ 

a In r 



(14, 



rflp " oi In r ' uk 

with fK Kcplcr velocity. As a simplification wc take the 
rate of change of the viscosity the same at all heights z and 
replace the logarithmic derivative of /x by the logarithmic 
derivative of its integral / 

f = 2[\dz (15) 

J Zn 



d In /i d\nf 



d\ni 



film 



(integration over the height of the corona on both sides 
of the disk). / results from the computed vertical distri- 
bution of the pressure according to Eq. (10). With this 
replacement Eqs. (10), (14) and (15) yield for the coronal 
accretion rate (positive if inward) 



Mc = 2 2nrpvrdz = 37r/(l + 2—-^) 
Jzo a m r 

= —2TTrT,Vr 



(16) 



the third equality defining the mean radial velocity Vr, 
with the surface density S = 2 J^^ pdz (both sides of disk). 



2.3.2. The divergence of mass and energy flow 

We derive the modifications of the sidewise advection 
terms in the equations for conservation of mass and en- 
ergy again as mean values over z. The radial divergence is 

written as 



1 d . ^_ . 2^_ 

- — [rY,Vr) = i^VrVM 

r ar^ ' r 



(17) 



such that r\M measures the radial divergence in terms of 
its one-zone model value. Thus follows 



1 Idln(rStJ^I) ldln|Mc 



rflnr 



2 dlnr 



(18) 



With Eq. (16) the mass modification factor can also be 
written in the form 



1 din 1/(1 + 2^ 



dlnr 



(19) 



The divergence of the energy advection involves the 
derivative of • u{r,rj) with u the specific energy (see 
below). 



ldl n|Mc(r,T7) ■M(r,T?)| 
rflnr 



(20) 



The potential, kinetic, and thermal specific energies 
all scale approximately as 1/r (see Liu et al. 2002, Fig. 4). 
The factor then is 



rjE = r]M + 0.5. 



2.3.3. The specific energy 



(21) 



(13) In the specific energy u the kinetic energy is ^ = (w^ 



+ ^7p)l^- Since v^^, the dominant term, is r dependent, 
and it has to be taken into account as vj^ = GMr"^ jir^ + 
2;2)3/2 ^jj^ ^j^g one-zone model it was taken as constant 
and dropped out of the equations). 

2.4. Boundary conditions 

The boundary conditions are the same as those in the 
one-zone model. 

As the lower boundary of our calculations at z = zq, 
we take the level where electron temperature and ion tem- 
perature are already about the same and have the value 

— Ti ^ 10^'^K. Below an analytic solution suffices to 
cover the small remaining extent before the disk chromo- 
sphere is reaching it and yields a relation between pressure 
and downward heat flux (for a more detailed discussion see 
Meyer et al. 2000). 

At the upper boundary with no pressure confinement 
at infinity we require sound transition at some height z = 
z\ (free boundary), Vz — V&. Further with no influx of heat 
from infinity there, we require Fc = neglecting a small 
remaining outward heat fiow.. 

This constitutes 5 boundary conditions for the 5 ordi- 
nary differential equations in z. 
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3. A consistent solution at all distances r 

The equations described are valid for accretion in differ- 
ent astrophysical objects, stellar and galactic black hole 
masses. The solutions scale with Eddington accretion rate 
and Schwarzschild radius. In the following we present the 
solutions for an accretion flow around a QMq black hole 
(the Eddington accretion rate then is 8.4 • 10^®g/s, the 
Schwarzschild radius 1.8 • lO^cm). We take the viscosity 
parameter a=0.3 (for the dependence of the coronal struc- 
ture on the viscosity see Meyer-Hofmeister & Meyer 2001). 

3.1. How to find consistent solutions 

The main goal of the new approach is to find out how the 
corona is stretched out above the disk in radial extent. The 
one-zone model gives a description of the coronal structure 
only for the one inner zone where the thin disk underneath 
becomes truncated. This is a reasonable first approach 
when the evaporation efficiency is steeply increasing with 
decreasing r and evaporation is the dominant mass feed- 
ing cause there. (Mass flow into the corona from farther 
outward was neglected and simple estimates for the flow 
of mass, energy and angular momentum through this one 
zone were made). 

For the determination of the coronal mass flow in a ra- 
dially extended region the difference of inflow and outflow 
in a succession of radial zones is now properly accounted 
for. This means that the coronal mass flow at a certain ra- 
dius depends on the coronal structure in the neighboring 
zones and must be determined iteratively. 

We first compute the coronal structure for successive 
distances and a series of values for the modification terms 
to construct a grid of curves /(r, 77M, ??£;)■ The derivative 
of / only appears in Vr and in the modification terms. We 
therefore take these terms together as 



Vm = Vm • (1 + 2 
Ve 77b • (1 + 2 



dln/(r,?/) . 

d\nr ■ 
dln/(r, 7?) . 

dlnr 



(22) 



We show some of these results in Fig. 1. A value 
f/M — means a net zero outfiow. In this case we found 
that the structure approximately is independent of the 
value oifjE- The value of the viscosity integral / decreases 
with increasing net flow. Higher values of / correspond to 
higher values of the pressure at the bottom of the corona. 
The value f/M = 1 corresponds to the net outflow taken 
in the one-zone model (aside from the now included term 
u^). From Fig. 1 we see that advection of energy has a 
different effect on the coronal structure at different dis- 
tances r. This depends on whether the frictional energy 
release preferentially goes into radiation or into advective 
loss (compare Fig. 5 Meyer et al. 2000). 



What we finally need is a series of coronal structures 
at each distance determined for given 77 values (77 indicates 
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8.0 



8.5 



9.0 



9.5 



log r 



10.0 



Fig. 1. Values of the viscosity integral /(r, t^m, 77^) for dif- 
ferent parameters, basic curves for the evaluation of the 
coronal flow, solid lines: (?7m, 77£;)=(0,0), (0.5, 0.5), (1,1); 
long dashes: (77M, ■^£;)=(0.5,1),(1,1.5) For (0,0.2) the line 
is approximately the same as for (0,0). Thick line: exam- 
ple of a consistent curve log /, including the mass, energy 
and angular momentum flow consistently (in our compu- 
tations we investigate the accretion flow around a 6Mq 
black hole). 

the dependence on both 77 jv/ and 77^) so that the curve log 
/(r, 77) has the slope s = din f /dlnr that flts to these 
values taken for 77. From Eqs.(19), (21), and (22) and the 
deflnition of s we obtain 



1 , 9 ds 



), 77B = (77M + 0.5)(l + 2s).(23) 



This implies an ordinary differential equation of second 
order for f{r). A unique solution then requires two bound- 
ary conditions, e.g. at an inner and an outer edge. The so- 
lutions thus form a two-dimensional manifold. Different 
approaches might be possible. An iterative procedure 
without enforcing appropriate boundary conditions, start- 
ing from one of our curves in Fig. 1, determining the 
derivatives of log / and the 77 values and then the ap- 
propriate new log / curves, seems not to converge to a 
consistent solution. This is probably due to the large num- 
ber of possible solutions corresponding to the diversity of 
possible boundary conditions. 

We use the following procedure. We start at a chosen 
distance r = rtr where the thin disk is truncated. There 
the boundary condition is s(7'tr) = which means that all 
angular momentum carried inward into the inner disk free 
region is returned by friction (see discussion in Meyer et 
al. 2000). Now we construct a consistent log / curve step 
after step. Step 1: we choose an initial value of log /(r, fj) 
(from our computed grid of curves); from Eq. (23) we see 
which value has to be chosen for -7^^ to give the values 

a In r ^ 

fjM and fjE that belong to the chosen initial value of / 
(the derivative is approximated by a difference quotient 
(s(logrtr + ^logr) - s{logrti)/ Alogr, we used Alogr=0.2). 
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This determines the slope s with which to proceed to the 
next log / value at the next radius. Step 2: the slope de- 
termined in step 1 yields the next value for log /; again 
we determine the further slope so that relation (23) is ful- 
filled. Each step thus determines the further outward slope 
s of the log / curve so that the ry values and /(r, fj) are 
consistent. 

For a given truncation radius and for each initial value 
of log / we get with this procedure one consistent log 
/(r, 77) curve. A second boundary condition at an outer 
boundary of the interval for which we determine the 
coronal flow allows to discriminate between the various 
curves that belong to one truncation radius. The slope at 
the outer distance is a possible boundary condition. For 
s = 1/2 the radial diffusive velocity Vr becomes zero, no 
mass flow inward or outward. Then closer to the black hole 
the mass flows inward, at larger distance the mass flows 
outward. In different disks the outer boundary condition 
might be different. In the case of X-ray binaries it seems 
plausible that tidal forces at the disk boundary prevent 
outer coronal mass outflow, that is din f /dlnr becomes 
-1/2 there. 

3.2. Computational results 

Our method yields the coronal structure solutions for 
given truncation radii. This gives the mass flow rate in 
the corona at the truncation radius where no mass flows 
in the cool disk anymore. This mass flow together with the 
wind loss from the corona added up over the radial extent 
of the corona is equal to the mass flow in the far out cool 
disk. This establishes the dependence of the disk trunca- 
tion radius on the mass flow rate in the outer thin disk. In 
Fig. 1 we show one log / curve starting at log rtr=8.5 In 
Fig. 2 we show three characteristic curves starting at log 
rtr=8.5 as well as the characteristic run of log / curves for 
other truncation radii. Generally the lower the initial of 
log / the steeper is the corresponding log / curve. One also 
sees that a small variation in the initial value of / yields 
a large variation of the radius where Vr = (marked by a 
diamond in Fig. 2). This indicates a somewhat weak de- 
pendence of the inner coronal mass flow rate on the outer 
boundary condition. 



Comparing the / curves of the solutions with the field 
which our basic grid curves cover in the log /-log r dia- 
gram we see that the solution /- curves can easily lead 
into areas outside of the grid, either to higher or lower 
values of log /. If they stay within the range of computed 
basic curves necessarily the slope becomes less than -1/2. 
No consistent solutions have been constructed which have 
Vr=0 at a very large distance. The too steeply decreasing 
curves correspond to the situation that the coronal flow 
is inward only near the inner boundary but turns into an 
outflow already a short distance away from it. To continue 
these curves correctly we need the coronal structures for 
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8.0 8.5 9.0 9.5 10.0 

log r 

Fig. 2. Consistent solutions for log / for truncation radii 
log r=8.1 to 9.3. Diamonds on two curves mark where 
Vr — Q. The coronal flow at the truncation radius is pro- 
portional to /. Dashed parts indicate uncertain extrap- 
olation outside of the computed basic grid curves. The 
dotted curve is identical with the uppermost curve in Fig. 
1 (accretion flow around a 6Mq black hole) 



high ffM values, which means a high net radial mass out- 
flow. The flatter curves also can leave the fleld of standard 
solutions, that is fjAi < 0. The corona would then have a 
negative radial net outflow, which means less mass radi- 
ally leaves the zone at the inner radius than enters into it 
at the outer radius. 

In Fig. 3 we show the coronal mass flow for the case 
of a thin disk with a mass flow rate of 2.4 • lO^^g/s from 
outside. This disk becomes truncated at log rtr=8.5. We 
show how the coronal flow increases inward (together with 
the mass and energy flow modification values and fjE oi 
this consistent solution). The evaporation rate measured 
by A/evap = ^irr^rriQ (toq vertical mass flow rate per unit 
surface at the bottom of the corona, in units of g/cm^s) 
gives the amount of gas locally evaporated into the corona. 
The earlier one-zone model only considered evaporation in 
the innermost region. Here our solution shows a radially 
quite extended evaporation region. In the lower panel of 
Fig. 3 we show how the mass flow in the thin disk decreases 
inward. Gas evaporates into the hot coronal flow, from 
which a signiflcant part, « 50% is lost in a wind. 

In Fig. 4 we show the results from a series of compu- 
tations with successive truncation radii. We chose solu- 
tions where ^^=0 occurs relatively far outward and give 
the coronal mass flow rates at the truncation radii. In 
addition we show the virial and ion and electron temper- 
atures of the corona at height z — r at the truncation 
radius. Note the small difference between Ti and T^. The 
coronal mass flow rates are higher than the ones derived 
with the one-zone model. This is in part due to inclusion 
of the value of v'i in the specific energy. 
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Fig. 3. Upper panel: solid lines: Mc coronal mass flow rate 
above the thin disk around a &Mq black hole, truncated at 
log rti— 8.5, mass and energy flow modification values t)m 
and f)E- Mevap evaporation rate (see text). Dashed lines: 
Mf, and Movap for disk truncation at log rtr=9-l- Lower 
panel: mass flow rate in the thin disk (gray area) and in 
the corona (dotted area) as function of distance r from 
the black hole, area remaining above indicates the rate of 
gas loss in the wind from the corona. 



4. Discussion and conclusions 

4.1. The nature of the solutions - wind loss 

Having developed the procedure for solving the differen- 
tial equations described above we see that there are quite 
different regimes of solutions. Up to now we have studied 
solutions with a moderate advective mass and energy net 
outflow. These solutions seem adequate to describe the 
coronal flow above a thin accretion disk e.g. in X-ray bi- 
naries. For the investigation of AGN disks another outer 
boundary condition might be taken, but this would not af- 
fect the solutions essentially. It will be interesting also to 
study other regimes of solutions i.e. solutions leading into 
the area above our network of curves log f{r,fiM,flE) in 
the log /-log r diagram. Test computations showed that in 
these solutions an increasing part of the coronal flow does 
not proceed inward but is lost in a wind. In our equa- 
tions flaring vertical columns arc considered. If the wind 
loss becomes important the structure becomes dependent 
on the assumed geometry. This geometry then has to be 
considered carefully to have a consistent picture how the 
expanding wind structure is arranged radially. This is an 
interesting topic in view of the question how much gas ac- 
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Fig. 4. Thick line: series of solutions for the coronal mass 
flow rate Mc, for disk truncations at log r from 8.5 to 9.7; 
thin lines: ion and electron temperature at coronal height 
z — r together with the virial temperature T = 



tually enters the disk free inner region and flows towards 
the black hole. 



4.2. Disk truncation closer to the black hole? 

If we imagine that the mass flow rate in the thin disk can 
be arbitrarily high we should encounter a situation where 
the disk is not yet truncated at e.g. log r=:8.5. If evapora- 
tion continues a truncation farther inward is possible. We 
have not yet found a consistent solution for the coronal 
flow corresponding to a disk truncation farther in. As can 
be seen from our solutions starting at larger distances we 
would have to start with higher initial values of log /, val- 
ues for which we have no log / curves. Logically, a coronal 
structure for condensation of gas instead of evaporation 
should then be appropriate. Such condensation already 
appears in the algebraic treatment of coronal structure 
by Rozanska and Czerny (2000a). Thus in our approach 
which consistently derives advective terms from the ra- 
dial gradients condensation should also be expected. In 
the case of a net mass inflow by sideways advection part 
of the gas flows down into the cool disk while another part 
flows upward and leaves as a wind. 

4.3. Compton effect 

The computed coronal structure is valid for a mass flow 
rate in the thin disk low enough not to affect the corona by 
Compton cooling. Results in the work of Liu et al. (2001, 
Fig. 5) have shown that for higher mass flow rates the 
Compton cooling leads to a steep decrease of the evapo- 
ration rate inwards. For our investigation this means that 
then the basic grid of log / curves (Fig. 1) would also have 
a more pronounced decrease towards smaller r. This is the 
region where we have not yet constructed consistent solu- 



8 



E. Meyer-Hofmeister and F. Meyer: The formation of the coronal flow/ADAF 



tions. Having solved the problem of solutions for smaller 
r it is meaningful to discuss the effect of Compton cooling 
on consistent coronal flow solutions. 

4.4. Spectral transitions in X-ray binaries 

The two examples for coronal flow corresponding to disk 
truncation at log rtr equal to 8.5 and 9.1 (see Fig. 3) show 
the difference in the mass flows in disk and corona. For a 
mass flow of about 9 • 10^''g/s the thin disk becomes trun- 
cated at log rtr =9.1. The hard spectrum is determined by 
this amount of mass flow in the corona, farther inward in 
the ADAF. If the mass flow in the thin disk is 2.4 • 10^*g/s 
the disk is truncated at log rtr= 8.5 (compare Fig. 3 for 
the amount of wind loss), the spectrum correspondingly 
changed. A factor of about 2.5 in the mass flow rate might 
occur during the increase or decrease to/from an outburst 
in an X-ray transient (even larger if wind loss is included) . 
This leads to the change of the spectrum. For very high 
rates the disk is not truncated. For this case we can not 
yet present a solution. 

4.5. Concluding remarks 

The newly developed method for the construction of ra- 
dially extended coronal flows allows a better description 
for the change from accretion via a thin disk to a disk free 
pure coronal flow (ADAF). We have presented results for 
disk truncation at various distances r from a stellar black 
hole of 6Mq applicable in X-ray binaries. The method also 
allows to study other regimes of solutions and might shed 
light on wind loss under those conditions. Particularly im- 
portant will be to investigate whether the corona can par- 
tially condense into the thin disk again or keeps flowing 
at a significant rate. The solutions scale with Eddington 
accretion rate and Schwarzschild radius, are therefore also 
relevant for the accretion flow in AGN. 
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